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Abstract — This paper considers the problem of channel sensing 
in cognitive radios. The system model considered is a set of N 
parallel (dis-similar) channels, where each channel at any given 
time is either available or occupied by a legitimate user. The 
cognitive radio is permitted to sense channels to determine each 
of their states as available or occupied. The end goal of this 
paper is to select the best L channels to sense at any given 
time. Using a convex relaxation approach, this paper formulates 
and approximately solves this optimal selection problem. Finally, 
the solution obtained to the relaxed optimization problem is 
translated into a practical algorithm. 

□ 

I. Introduction 

As the number and types of wireless (multimedia) applica- 
tions increase, so do the stringent requirements they impose on 
the wireless medium. Thus, it is essential that we determine 
efficient means of utilizing limited spectral resources avail- 
able to us. Currently, bandwidth resources are divided into 
frequency bands and allocated to different users exclusively 
in order to insure the quality of service (QoS) of multiple 
wireless systems, and the FCC's frequency allocation chart 
[1] shows that almost all frequency bands are currently divided 
and allocated to different groups for varying purposes. Also, 
according to recent surveys [2] and [3], most of this allocated 
radio frequency spectrum is vastly under-utilized. Cognitive 
radios are emerging as promising solutions to enable better 
utilization of spectrum especially in bands that are currently 
underutilized [4]. The classical example of a cognitive radio is 
one that employs "interweave" cognition [4]. These interweave 
cognitive radios are permitted to occupy a channel (frequency 
band) only when it is not occupied by a user licensed to 
use that band. If the presence of other radios can be sensed 
accurately and quickly, then such a policy can help ensure that 
cognitive users cause little to no interference to the licensed 
radios in the system. A majority of existing literature on 
cognitive radios focuses on such interweaved radios. For an 
analysis of other classes of cognitive radios, see [5], [6] and 
[7]. One of the main issues under study in the interweaved 
cognitive radio domain is the so-called "sensing problem", 
where we desire to determine, as accurately and efficiently 

'This work is supported by a grant from Samsung Advanced Institute of 
Technology. 



as possible, if a given channel is occupied at any given time 
[8], [9], and [11]. For example, [8] describes a simple energy 
detection scheme for additive white Gaussian noise channel. 
In [9], the performance of energy detection schemes in a 
multipath environment is analyzed, and in [11], the impact of 
additional side information is considered in determining the 
performance of cognitive sensing. Overall, channel sensing is 
one of the better established fields of research on cognitive 
communication. In this paper, our goal is significantly different 
from that of channel-sensing literature. Given a fairly accurate 
sensing algorithm, we desire to determine which channels 
should be sensed when. In addition, we desire to perform a 
resource-allocation problem across multiple channels which 
may or may not be available to the cognitive radio. Overall, 
we ask the question "Given that there are multiple dissimilar 
channels available for you to sense, which channels should 
you sense and, if they are available, what rate/power should 
you assign to them?" 

The dissimilarity between different channels arises from 
various factors. The properties of the propagation environment 
depend on frequency and thus can be significantly different 
from channel to channel. Some channels may suffer from 
"extraneous interference" from non-legitimate sources (such 
as in the industrial, scientific and medical (ISM) bands) that 
reduce the channel quality. Thus, just as any other multiband 
radio, the cognitive radio must allocate resources across dif- 
ferent bands it uses while simultaneously determining which 
ones it is permitted to exploit. Note that, in isolation, the 
problem of channel selection for cognitive radios [12], [13] is 
well studied. Also, by itself, the resource allocation problem 
for muti-band radios is also well-understood [14]. However, 
bringing the two together is both important and challenging as 
they are tightly coupled in the context of interweaved cognitive 
radios. A simple explanation of this strong interdependence 
between sensing (channel selection) and resource allocation is 
as follows: Let us say that the system is such that "noisy" 
(poor) channels are less frequently used by licensed users 
than "clear" (good) channels. If the sensing mechanism were 
to choose to sense the infrequently-used channels, it will 
present the cognitive radio with available channels that are all 
"poor" resulting in a low rate. On the flip side, if the resource 
allocation mechanism were to assign high rates to the "good" 
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Fig. 1. Channel Model 



channels, the sensing mechanism may find that they are not 
available for use and then again sustain a very poor rate. Thus, 
designing channel selection and allocation jointly is essential 
for cognitive radios. Note that this paper's focus is on the 
fundamental limits of joint selection and resource allocation 
in cognitive networks to provide a benchmark on performance. 
Thus, aspects such as sensing error, delay, device and network 
non-linearities etc. are not incorporated into the analysis. 

The rest of this is organized as follows. The next section 
details the system model and notations used in the paper. 
In Section III, we find the fundamental limit of the given 
system model. In Section IV, we propose an algorithm for 
joint channel selection and power allocation, and we conclude 
with Section IV. 

II. System Model and Problem Statement 

The channel model is shown in Fig. 1. We consider N 
parallel legitimate channels with equal bandwidth. In each time 
slot, a channel n, 1 < n < N, is occupied by a legitimate user 
with probability q n . There are one cognitive transmitter and 
one cognitive receiver. The cognitive transmitter is allowed to 
transmit over channel n, if it is not occupied by any licensed 
user. In legitimate channel n, cognitive radio's channel is: 

where Z n is additive Gaussian noise of variance cx^. Note that 
this noise variance can be different from channel to channel, as 
it represents the fading state of that particular channel. Before 
the start of cognitive radio's transmission using the legitimate 
channels, the cognitive transmitter should know whether they 
are occupied by the licensed users or not. Thus, at the start of 
every time slot, the cognitive transmitter is allowed to sense 
a subset of channels, and is allowed to exploit those channels 
that are unoccupied; in this paper, we assume that the sensing 
is performed perfectly. Also, the cognitive transmitter is not 
allowed to transmit using the channel which is not sensed in 
order to guarantee the transmission of the licensed users. As N 
is assumed to be large, it is impractical to allow the cognitive 
radio the ability to sense all of them at the start of every slot. 
Instead, we require it to cleverly choose a subset of bands on 
which to focus its efforts. The capacity of the cognitive radio 
depends on which channels to sense from N parallel channels, 



and power allocation among the available parallel channels. 
Average total transmission power of cognitive transmitter is 
constrained to P. 

First, define the I n (t) and lE.n(t) to be the indicator 
function for selected channel to be sensed and an indicator 
function for the unoccupied channel respectively, i.e., 
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and 



if channel n is not to be sensed 
if channel n is to be sensed 

(1) 

if channel n is occupied 

1 if channel n is unoccupied 
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Denote the time average capacity of the cognitive radio with 
the selection of the sensing channel I n (t) and power allocation 
P n (t) in one time block as C (I n (t), P n (t)). Then, 
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where T is the number of time slot in each time block. 

In our model, we assume two constraints on the cognitive 
radio: 

a. An average power constraint on the cognitive transmitter of 
P, 

b. The number of channels that can be sensed by the cognitive 
radio at any given time is L < N. 

Note that if the cognitive radio could sense all channels, L = 
N, this problem has a fairly trivial solution. At the start of each 
time slot, the cognitive radio would determine all available 
channels and waterfill its power over them [10]. 

If the number of channels that cognitive radio can sense is 
less than N, i.e. L < N, the resulting optimization problem 
can be stated as follows: 



max C (/„(*), P„(t)) 
P„(t),i n (t) 



such that 
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Pn(t) > 0, 
/n(t)G{0,l}, 

i B>n (i)e{o,i}- 



(4a) 

(4b) 
(4c) 

(4d) 



The optimization problem given by (4) determines the 
maximum empirical average rate achieved by the cognitive 
radio given constraints on the system. Note that it is an integer 
programming (IP) due to the constraints in ( |4d| , and multi- 
dimensional due to its dependence on time t. 

The next section studies the optimization problem given by 
(4) in an ergodic policy setting. 



III. Optimal Power Allocation and Selection of 
Sensing Channel 

As a first step, we assume that our policy is ergodic and 
"static" , i.e., that our sensing and power allocation policies are 
only functions of the channel statistics and do not evolve with 
time. This results in the following (simplified) optimization 
problem: 



max 



such that 



N 

E 



Inq n , A . Pn 

—— lOff 1 + -7T 



J2 ^VnPn < P, 
n=l 

N 

E 7 » ^ x < 



(5a) 



(5b) 



(5c) 
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and 



(5d) 



Pn>0, 
In G {0, 1}. 

Since P n = where J„ = 0, constraints pb] ) can further be 
relaxed to 



N 



InPn < P- 



(5e) 



Denoting the optimal selection of channels to be sensed and 
power allocation for channel n as I* and P* respectively, the 
optimum solution for (5) is given by the following theorem: 

Theorem 1: The joint channel selection & rate allocation 
problem (characterized by the optimization problem in (5) is 
maximized when: 
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where 
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N 

E 7 « = L ' 



and \w] + is maximum value of and w. 

Proof: Note that ( f5a] ) is a concave function over P n for a 

particular choice of /„. The following Lagrangian describes 



the optimization of (5ai with respect to P n for a given /„: 
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Taking the derivative of ([6]) and setting it to zero, we get 
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Fig. 2. power allocation with given (I\ = 1, 1% = 1, ^3 = 0, Iff = 1) 
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From dSeb we obtain, 
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[A - cr^] + I n q n = P. 



(9) 



Fig. 2. provides a graphical representation for the power 
allocation strategy in dHl. Note that it is similar to the water- 
filling solution, with the main difference that the each channel 
has different width, q n I n . We refer to the policy in (|9|l as 
modified water-filling throughout this paper. Given that we 
understand the structure of the power allocation policy that 



optimizes (5a 1, we now desire to determine I*. Note again that 
the optimization problem in (5) with respect to /„ is an IP. 
It can be found by an exhaustive search, but computationally 
very hard to solve. Moreover, the power allocation strategy in 
[S] specifically, the water-level A is tightly coupled with the 
choice of I n . In the next section, we present an algorithmic 
framework that approximates /* (and thus the water-level A) 
using low-complexity iterative techniques. 

IV. Joint Selection and Power Control 

A typical integer program is non-polynomial in complexity. 
Although multiple techniques exist for obtaining approximate 
solutions to such a program (such as branch and bound [15], 
relaxation), such techniques apply to any integer program and 
do not take the structure of the problem into consideration. 
Our focus is on developing an algorithm customized to this 
problem setting. We perform this in two steps, which we call 
"coarse" and "fine" optimization. The coarse optimization step 
determines a set of L channels to be utilized by the cognitive 
radio. It gives us the lowest possible waterlevel, A m ;„. The fine 



optimization step uses \ m i n which we obtained from coarse 
optimization to further optimize the choice of the L channels. 
First, we describe the coarse optimization step: 

Coarse Optimization: We iteratively find the channels to 
sense along with modified water-filling which incur the lowest 
water level. Let X m i n denote the lowest water level, and 1° and 
P£ indicate the selection of the channel and power allocation 
which result in \ m in- Detailed procedures to find \ m in, 1%, 
and P° is described in the following four steps. 

Step I: Start with L initial channels. We can choose L 
channels with the largest q n as initial channels, for example. 



, 1 if q n is among L largests 
ln '°~^ otherwise (1U) 



then the coarse optimization is optimal. 

Proof: Define S* to be the set of channels from optimal 

selection. 



From definition, 
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.7 = 1 



(11) 
(12) 



Step II: Perform the modified water-filling with J n j > 1, 
such that 



JY 



E \ X i ~ ^j-iln = P- 



(13) 



n=l 



Step III: Calculate q n (Xj — a 2 ), and select the largest L 
channels. 
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n ,j — \ Qn(Xj — <7„) is among L largests 

otherwise 



Let's assume that there exist at least one legitimate channel 
with noise variance higher than the lowest water level which 
is included in the optimal channel selection. 

3n G S* , > A min . 

Define S' — S* U S c , and allow the number of channel to 
sense to be M = \S* U S c \, which are strictly larger than L. 
Then, 

max flog (l + > m - E | l0g ( 1 + § 

n£S' V ™ 7 riGS* V ™ 

(20) 

Note that 5 C C S", and cr^ > \ min for all rt ^ S c . Modified 
water-filling of M channels in S' will lead to P n =0 for all 
n ^ S c . Thus, 

max ^ ^ log ( 1 + | = max Gr log I 1 1 



Combine the above result with (1201, we obtain 



n£S c 



(14) 



S,- = {ne[l,iV]|J n)J - = l} (15) max^ |log(l + ^) >max^ |l g(l + ^) 



Step IV: If Sj — Sj—i, terminate the iteration, and set the 
power allocation and channel selection values. 



Xmin X n 

r = i 



(16) 
(17) 
(18) 



Otherwise, j = j + 1 and repeat from step II. 

The coarse optimization is performed for two reasons. One 
is that the performance of coarse optimization is very close 
to the optimum. This will be shown from the simulation 
result in the next section. Here, the optimality of the coarse 
optimization in one special case will be stated and proven. 

Lemma 1: Define S c to be the set of the channels which 
are selected from coarse optimization; 

S c = {ne[l,N]\I* = l}. 

If the noise variances of all the channels which are not selected 
in the coarse optimization are greater than the lowest water 
level X min , i.e., 

nun : 



(21) 

Above result contradict ( p~9] >, unless S c is the optimal. This 
concludes the proof. 

The other reason for performing the coarse optimization 
is that it provides the essential information, X m in, which is 
necessary for further fine optimization. Upon the following 
assumption, fine optimization is optimal. 

Conjecture 1: If the noise variance tr^ is greater than the 
water level in the coarse optimization (er^ > A m j„), then the 
channel n is not likely to be sensed in the optimal strategy, or 
even if it is included it will not increase the capacity much; 
Intuition:The following gives the intuition for the above 
conjecture. Define to be the set of channels with noise 
variance greater then or equal to the X m i n and S~ to be the 
set of channels with noise variance less then X m i n but not 
included in S c ; 



S+ = {ne[l,N]\a 2 n >X min }, 
S- ={ne [l,N]\o* < X min ,n 



S c }. 



We have the set of channel S c which incur the lowest 
waterlevel. Lemma([T]) shows that average capacity cannot 
increase by exchanging elements in S + with elements in 



S c . Thus, for elements in S + to be included in S*, optimal 
channel selection, elements in S~ should be included also. 
By exchanging elements in S~ with elements in S c the 
waterlevel rises up. Elements in S + can only be in optimal 
selection if exchanging S + with elements in S c lower the 
waterlevel which increased due to the inclusion of channels in 
S~ effectively. Intuition is that channels in S c are the channels 
which can lower the waterlevel effectively already. Thus, effect 
of lowering the waterlevel with channels in S + will not affect 
much in increasing the average capacity. The validity of this 
conjecture is shown from the numerical analysis. 

Fine Optimization: From lemma([T]), if the number of chan- 
nels that is selected to sense from the coarse optimization 
is less than L, it is optimal, and no further optimization is 
needed. Otherwise, further optimization will be required. From 
Conjecture^, we reconstruct the problem, so that we can 
optimize the selection of the channel over the channels with 
noise variance smaller than or equal to lambda m in only, we 
rearrange the useful channels by indexing from 1 to M, where 
M is the number of channels that has noise variance smaller 

than \ m in\ 

M=\S c US-\ (22) 
a 2 n -X mm /leO VnG [1,M]. (23) 
Then, the optimization problem can be rewritten as follows; 
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such that 
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f A - °l] + ^Qn = (A - al) I n q n = P, (28) dl n 



(24) 
(25) 

(26) 

(27) 



A I 



In e {0, 1}, 



(29) 

(30) 
(31) 



where (a) and (b) result from constraining A > X m %n- Then 
the optimal channel selection and power allocation can be 
determined by using the following theorem. 

Theorem 2: 



A 



£n=l 1nIn<J n + P 

EM j 



I n = if A > ale l ~^ 



Proof: Relax the constraint on /„, such that the can take 
the value in the region [0, 1]. Construct the objective function 
C (A, / (I n )) such that it is concave over the region of I n 
and A, and /(0) = and f(l) — 1. Consider the function 
/ (/„) = /*, then /(0) = 0, /(l) = 1, and 



C (A, /(/„)) 



M a T k \ 



n=l 



log 



(32) 



Concavity of C (A, / (I n )) can be found as follows; 



a 2 e(A,/(/„)) 
oil 

3 2 C(A,/(/„)) 
d\dl n 

q n k (k - 



3 2 C(AJ(/„)) 
dl n d\ 

a 2 c(A,/(/„)) 
a 2 A 

n/*-2i„„ a 



(33) 



log^ 



q n kl k n ^loge 



"loge 
E^li^-^loge 



(34) 



Since the matrix is symmetric, if the determinant, (1, 1) and 
(2, 2) components of the matrix take the negative values, the 
matrix is negative semi-definite. Thus, 



k(k - 1) log ■ 



k(k- 1) < 
f k 2 log e < 



(35) 
(36) 



are the condition for C (A, / (I n )) to be a concave function. 
We can find k such that the condition can be satisfied; 

/ bg A_ v 
k = min I : , ^ ) . (37) 



log 



log e 



Now that we verified the concavity of the objective function, 
we can construct the according Lagrangian multiplier: 

Ef =1 ^iog(4) 

C = ~Mo (Ef=i fei?(A - °l) - P) - Mi (Ef=i - £) 

+ E„=1 M2,i-fi - E„ = i M3,i(Xi - 1) + /m(A - A mm ). 

(38) 

Solving the optimization, 
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+ M4 = 



\q n I k n (\-a 2 n )~P \ =0, 
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H2,iln — 0, 



(39) 
(40) 

(41) 

(42) 
(43) 



Ms.iUn - 1) = 0, 



(44) 



fj,4,{\ — A m i„) — 0, 



(45) 



where /xo, Mi> M2,i> A*3,i> an d M4 are non-negative values. From 
the condition d40l, 



A = 



MO E„=l ?n J n - M4 
„ rfcl v^Af 



(46) 



Note that^^ =1 g„J^loge, J2 n =i Qnl%, and M4 are the non- 
negative value. Thus, jUq should be positive number in order 



for A to be positive. Then, from the condition (41 

M 
n=l 

Thus, 
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and from rearranging the equation (|46 
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From the condition d39|l, we find that 
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(49) 



qn log - M0<7n(A - <) - Ml ) fc 



M3,i — M2,' 



(50) 



If 7 n is not either or 1, from conditions ( |43] l and ( |44] >, 
/i2,i and /i3 i becomes 0, which will result in making I n 
infinite number. Thus, I n takes either 0,1 value, which gives 
the desirable solution, such that the optimization of the re- 
laxed condition coincides with the condition for the original 
problem, and 



Eni 1 Qnln<% + P 



A = 



EM j 



(51) 



We set /X4 to be zero, then from (|4"9"|) and d50l), we obtain 



9n log - Ml ) fc 



M3,i - M2 

As a result, /„ can be 1 only if 



A > ale 1 ' 



(52) 



(53) 



This conclude the proof. With the TheoremQ, we can design 
iterative algorithm to find the optimal selection of channels to 
sensed iteratively. 

Step I: Set the channels from coarse optimization to be the 
initial channels. 



I, 



71,0 



1 i f n G S c 
otherwise 



(54) 



S = {ne[l,iV]|J n ,o = l} 

i = i 



Step II: Calculate the waterlevel Xj from (51 

EM j 2 



A, 



P 



EM j 



(55) 
(56) 



(57) 



Step III: Calculate Xj - of t e 
channels. 



(47) In,j = 



1 if 



A; 



cr„e 



otherwise 



, and select the largest L 

is among L largests 
(58) 



Sj = {n€[l,N]\I ntj = l} 



(59) 



Step IV: If Sj — Sj-i, terminate the iteration, and set the 
power allocation and channel selection values. A-^, X,{, and 
pf are the waterlevel, channel selection, and power allocation 
from fine optimization, then 



X f = A„ 
Pi 



(A, 



(60) 
(61) 
(62) 



Otherwise, j = j + 1 and repeat from step II. Following 
from theorem(??) this algorithm gives the optimal selection 
of the channels to be sensed and power allocation, with the 
assumption that the channels with noise variance greater than 
X m i n do not affect the optimization. 

V. Numerical Analysis 

In this section, we present numerical example of capacities 
for coarse and fine optimization along with optimal solution. In 
this example, dissimilarity between channels is implemented 
by adapting the multi-path fading, which will incur frequency 
selective channel. Also, occupation of the legitimate channel 
is modeled by having q n to be uniform in [0, 1] and identi- 
cally distributed. Sixteen legitimate channels are considered, 
N = 16, and cognitive radio is allowed to select and sense 
eight channels from all the legitimate channels, L = 8. Fig. 
3. compares the capacities of sub-optimal algorithms with 
the optimal one, where the performance of optimal channel 
selection is obtained from the exhaustive search. The graph 
shows that performance of the fine optimization meet with that 
of optimal one. Thus, it can be stated that the Conjecture ([T| in 
fine optimization is valid. Coarse optimization also performs 
optimally in the low SNR region. In the low SNR region, it 
is likely that <r^ > X m i„, Vn € [1, N],n ^ S c , because there 
are not much power to waterfill. From the Lemma ([TJ, coarse 
optimization is optimal in such case. It is also worthwhile to 
note that coarse optimization perform as well as the optimal 
one. 
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VI. Conclusion 

In this paper, fundamental limits of interweaved cognitive 
radio has been verified. In the case that there are large number 
of legitimate channels and only limited number of them can be 
sensed, the capacity has been analyzed. However, it requires 
exhaustive search over the combination of all the legitimate 
channels, which is not practical in terms of complexity. Thus, 
two steps of sub-optimal solutions have been developed. 
Coarse optimization is developed, and verified to be optimal 
in the low SNR cases, and further optimization is performed 
to ensure the performance in the high SNR region also. 
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